// $Id: mat.h 78 2012-07-20 21:30:18Z alan.watson.f $

#include <gsl/gsl_matrix.h>

#define mat_t			gsl_matrix

#define mat_alloc		gsl_matrix_alloc
#define mat_free 		gsl_matrix_free

#define mat_elt(m, iy, ix)      (*gsl_matrix_ptr((m), (iy), (ix)))
#define mat_ny(m)		(m)->size1
#define mat_nx(m)		(m)->size2

#define mat_fwrite		gsl_matrix_fwrite

#define mat_set                 gsl_matrix_set_all
#define mat_add                 gsl_matrix_add
#define mat_sub                 gsl_matrix_sub
#define mat_scale               gsl_matrix_scale
